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transmission capacity between different utilities to build probabilistic distributions for the spot-market prices 
and for electric trades in the system; and 

outputting probabilistic distributions for spot prices and trades at a plurality of time periods of a planning horizon 
and for a plurality of delivery points to assist in said supply management, 

5 

2. The computer process of claim 1 , further comprising the step of defining different scenarios wherein the maximum 
generating capacity of each utility changes under said different scenarios. 

3. The computer process of claim 1 or 2, further comprising the step of using the input weather forecast to create a 
10 probabilistic distribution of load forecasts for each utility of interest. 

4. The computer process of any preceding claim, wherein the input price functions are given implicitly. 

5. The computer process of claim 4, wherein a price function of a utility is approximated by a quadratic function which 
'5 is convex, positive and increasing. 

6. The computer process of any of claims 1 to 3 wherein the input price functions are given explicitly as a cost function 
of the generating units. 

20 7. The computer process of any preceding claim, wherein the probabilistic distributions for the spot-market prices 
and for electric trades in the system are created by sampling from electric load, price functions and cost distributions 
to create a sequence of deterministic models. 

8. The computer process of claim 7, further comprising the step of minimizing the cost of operating an electric utility 
25 system under each sampled set of price functions, transmission costs and electric loads. 

9. The computer process of claim 7 or 8 wherein the probabilistic distributions for spot prices and trades are deter- 
mined by a bucketing process comprising the steps of: 

30 creating ranges of spot prices; 

counting the number of spot prices that fall within each range and normalizing the counted numbers to obtain 
a probability for each range; 

and generating a histogram of the probabilities of each range as an approximation of the probability distribution 
of the marginal cost at each delivery point at a given time. 

35 

10. A system for use in an electricity network to manage power supply by forecasting spot-market prices of electric 
power and trading transactions at different delivery points, said system comprising; 

means for inputting maximum generating capacities of each utility competing within the market, price functions 
40 at which a utility is willing to sell its power at a given time period, a weather forecast, and a transmission cost 

and transmission capacity in different power lines of the electrical system; 

means for using probabilistic distributions of electric load and the price function of each utility and transmission 
iost and transmission capacity between different utilities to build probabilistic distributions for the spot -market 
prices and for electric trades in the system; and 
45 means for outputting probabilistic distributions for spot prices and trades at a plurality of time periods of a 

planning horizon and for a plurality of delivery points to assist in said supply management. 
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Note that function block 54 of Figure 5 calls the process of Figure 6 which performs the sampling. The sampling 
process is performed using a uniform random number generator. Uniform random number generators are available on 
most computers and under any programming language. For example, the C programming language provides the func- 
tion rand as part of its standard library header. See B. W. Kernighan and D. M. Ritchie, The C Programming Language, 

s Prentice Hall, 1988. The probability of the different scenarios, p/ is used in conjunction with the electric-load proba- 
bilities, n t k t to come up with a sample. The process of Figure 6 receives from the calling process (function block 54 of 
Figure 5) the value of t In function block 601, a random number x between 0 and 1 is sampled from the uniform 
distribution. Function block 602 defines a variable, y and sets its value to zero. The variable y is used to compute the 
cumulative joint distribution at any given /cand s. The value of the scenario counter, s, is set to 1. Function block 603 

10 sets the value of k to 1 . In decision block 604, the value of the cumulative probability, y, is compared to x. If x is less 
than or equal to y, then the functions ff and loads o* are returned to the catling routine (function block 54) to be used 
in the optimization problem. Otherwise, the value of k is incremented by 1 in function block 605. Decision block 606 
checks if /chas reached its maximum limit of K t If k<K p then the value of the cumulative probability, y, is increased by 
p^rc* -1 and the execution is sent back to decision block 604. Otherwise, the value of s is increased by 1 in function 

is block 608. Function block 609 updates the value of y by adding p^ 1 tc*" 1 to it. The execution then loops back to function 
block 603. Function block 610 is reached from decision block 604 when the value of x is less than the cumulative 
distribution at the point (k,s). 

System's Output 

20 

At termination of the simulation runs, we have the spot prices or marginal costs, \ e j t at each point /, t=1 ,...,/?, and 
at each time period of the planning horizon. Here, the superscript £ indicates the sample index. We also obtain the 
power flow, / , in the different transmission lines of the power network. Note that if yf . . , is greater than zero, then 
node / is transmitting power to node / under sample £ of our simulation. Of course, if y^ ^0, then y*. ^ ^0. . To get 
25 the probabilistic distribution of the spot-market prices, X € jt , we perform the following steps! 

Create ranges for X € jt . For example, we may want to have intervals of $5/MWH width. That is, we have the intervals 
(0-5], (5-10], (10-15] 

Count the number of \ € jt that falls within each range. We normalize these numbers to obtain a probability for each 
30 range. 

The resulting histogram is an approximation of the probability distribution of the marginal cost at delivery point / at 
time t. 

The previous process is called bucketing and is described in the flow chart of Figure 7. The process requires a 
35 number of buckets (i.e., ranges) which is defined by the user and the lower limit, u y -, of each range. We can assume 
that there is an extra range (after the last one defined by the user) that covers the range up to infinity. Function block 

71 sets the number of observations (so far) in each range or bucket to 0. The process then moves to function block 

72 in which the value of X e jt is compared against the values of u } to find the index y for which u£ )J jt <u y+1 . The value of 
the counter m ; - is incremented by 1 in function block 73 to account for sample / Function block 74 increases the value 

40 of € by 1 ; i.e., we are ready to process a new sample if any. Decision block 75 checks if I is less than L If this is the 
case, then we have more samples to study and the process moves back to function block 72. Otherwise, we have 
completed classifying into the appropriate buckets. Finally, function block divides the frequency, m y by L to obtain 
the relatiVe frequency for each range or the probability. 

We can follow similar steps to create probabilistic distributions for the amounts of power traded, / f , across the 

45 different power lines, (/,y). The resulting distributions can be helpful to the decision maker in deciding the amount of 
power to be generated and its price. For a power broker, these distributions can be used to sign buying and selling 
contracts with different suppliers and consumers. 

so Claims 

1. A computer process for use in an electricity network to manage power supply by forecasting spot-market prices 
of electric power and trading transactions at different delivery points, said method comprising the steps of: 

55 inputting maximum generating capacities of each utility competing within the market, price functions at which 

a utility is willing to sell its power at a given time period, a weather forecast, and a transmission cost and 
transmission capacity in different power lines of the electrical system; 

using probabilistic distributions of electric load and the price function of each utility and transmission cost and 
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to form a distribution of the spot prices and trades. Given that the number of combinations at each time period is SfXK p 
it may be hard to perform all of the previous calculations in a reasonable time. Instead, one can sample from the electric- 
load and cost distributions to create a sequence of deterministic models. The sampling process (i.e., simulation) is an 
approximation of the exhaustive one described above. The sampling process is repeated until a reasonable approxi- 
mation is reached of the distributions of trades and spot prices or until a maximum number, L, of samples is chosen. 
Consider the market equilibrium at time t under a sample of price functions f h h 1 , . . . , n, transmission costs hi , . . 

n, /=1 n, and loads d h a=1 n. We drop the indices t, s and k\o simplify notations. Due to the competition resulting 

from the bidding process, we expect the spot price to be the minimum cost at which the demand can be met while all 
suppliers satisfy their profit margins. In other words, one needs to minimize the cost of supplying power, 



ming f.U^gljc, jYl j . (2) 

The decision variable x,- represents the amount of power produced by supplier / while y f j represents the amount of 
power flowing from /to / Note that if two nodes are not connected directly with a transmission line, then the capacity 
of this line is set to zero. The minimization of (2) must be solved subject to Kirchhoff-law which requires that the total 
power flowing into a node must be equal to total power flowing out of this node. Mathematically, this condition is written 
as 



n n 

J^yj.i+Xi = ]Lyi.j + <1i' i = l n (3) 

where d t is assumed to be zero for nodes that do not have electric load and maximum generation, G b is zero for nodes 
that do not have any generation. Tosatisfythe transmission capacity, we require that O^y^u^-, ibi,...,n 1 /=i,...,n. Also, 
given the maximum generating capacity of utility /, we must satisfy the constraint 0<x-<G,-, One can enforce 

other constraints such as the maximum capacity of a node. 

Due to convexity, the minimization of (2) can be solved easily. See, for example, M. S. Bazaraa and C. M. Shetty, 
Nonlinear Programming: Theory and Algorithms, John Wiley & Sons, Inc., 1 979. If the price functions, /■, are quadratic, 
one can use a speedy quadratic programming approach as described by R. W. Cottle et al. and K. G. Murty, ibid. Note 
that by solving (2) one obtains the amount of power to be generated and transported from each supplier to each demand 
point. One also obtains the dual variables or the Lagrange multipliers associated with the constraints of equation (3). 
The dual variable associated with delivery point / represents the marginal cost of power at point /. That is, it is an 
estimate of the spot-market price at / under the given demand and price functions. 

We repeat the previous process by sampling a new set of price functions and electric loads then resolving the 
model. The sampling process continues until changes in the distribution of spot prices and the decision variables x,- ■ 
are within a pre-specified threshold. Another possibility is to repeat the process for a fixed number, L, of sampler, 

The flow chart of Figure 5 describes the sampling process. It starts in function block 51 in which the time index, t, 
is set to 1. The process then proceeds to function block 52 in which the sample counter, €, is set to one. Note that 
wheneveY we select a new sample, i.e., price functions, transmission costs, and transmission capacities, the value of 
£ is incremented by one (function block 56). In decision block 53, we check if we have reached the maximum number 
of samples, L If the answer is yes, i.e., £=L, then the process proceeds to decision block 56 in which we move to a 
new time period in the planning horizon. If €<L, then the process moves to function block 54 in which a new set of price 
functions (one for each utility) and electric loads at the different nodes of the network is sampled. This step, function 
block 54, calls the sampling procedure of Figure 6 which returns a set of price functions, electric loads, and transmission 
capacities and costs. The data are then used to construct the problem of minimizing (2) subject to the constraints of 
equation (3) and other system constraints. We do not discuss the solution process here. Any commercial package that 
supports quadratic programming could be used. As an output for the optimization, we obtain the amount of power 
flowing in the different lines of the network, / {nr and the marginal costs, n € jt , corresponding to equation (3) at the 
different nodes of the system. The process is tfien directed to function block 53 in which a new sample is generated. 
After reaching the maximum number of samples in decision block 53, the process moves to function block 56 in which 
the time counter, f, is incremented by 1. In decision block 57, if t<T, then the process goes back to function block 52; 
otherwise, the process terminates. At the end of the process of Figure 5, we have for each time period a set of L values 
for the spot-market prices at nodes t=1,...,n. We also have a set of possible power flows between the different nodes 
in the transmission system. 
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is then obtained by adding a linear term that takes into consideration the profit margin resulting in ff t . 

To reflect the importance of each scenario and to make the model more complete, a probability is associated with 
each scenario s. That is, a probability, p t s , is associated with each price function that reflects the likelihood of its oc- 
currence in the future. The value of p t s can be approximated using the frequency of each of the functions fj t ft 
from the past. Of course, the user of th^ tool may adjust these probabilities using his or her intuition if there is a need 
to do so. 

3. Using the weather forecast, create a load forecast for each utility ot interest. Usually, these tools require 
historical data for the load on the system of a utility and the corresponding weather conditions, i.e., maximum and 
minimum temperature during the day. Table. 2 shows an example of this data. 

Table 2: 



Load vs. Weather Data 


Date 


Max 


Min 


1 


2 


3 




5/25/95 


65 


43 


2101 


2090 


1045 




5/26/95 


68 


44 


2173 


2120 


1932 




5/27/95 


74 


53 


2318 


2292 


2056 





20 



25 



30 



35 



40 



45 



50 



The first column is the date, second column is the maximum temperature over the twenty-four hours of that day, and 
the third column is the minimum temperature during that day. The next twenty-four columns (only three are shown) 
provide the electric load in MWH at each hour of the day. By passing this data to a forecasting tool, it can create a 
model that correlates the maximum and minimum temperatures to the load. These tools are also capable of incorpo- 
rating the effect of time of day and day of year into the forecast. The historical data of any utility can be obtained in an 
electronic form from the power pool. In the future, the electric load of a delivery point, i.e., geographic region, will be 
posted on the World Wide Web. We do not discuss the forecasting process in this document. See, for example, R. H. 
Shumway, Applied Statistical Time Series Analysis, Prentice-Hall (1988), and A. C. Harvey, Forecasting, Structural 
Time Series Models and the Kalman Filter, Cambridge (1991 ). Note that the forecast is given usually for each hour of 
the coming week. The forecasting tool also provides a probabilistic distribution, i.e., mean and standard deviation, for 
the load at each hour. That is, using the weather forecast, one can obtain for each utility a probabilistic distribution for 
its load at each hour of the planning horizon. If the forecasting tool provides a continuous distribution for the load at 
each time period, the continuous distribution can be approximated by a discrete one. The number of discrete realizations 
at time f is assumed to be K t Each demand realization, , i=1 has a probability of rc t k that reflects the likelihood 
of its occurrence. 

The weather forecast is usually available on the World Wide Web and can be imported by our system easily to 
pass it on to the forecasting system. 

4. Transmission fi , in the different power lines of the electrical system. Our notation permits changing the 
transmission cost from one period to another. It also allows using multiple cost scenarios depending on the status of 
the electric system. We assume that for each transmission line (/,/) we know the maximum transmission capacity MWH 
which we denote by if f . If a power line is not available during a certain period, we can set its capacity, u^. ( , to 
zero. Note that we define a power line by its ends or nodes / and / To simplify the task of representing the electric 
network and its delivery and generating points, we assume that there are n nodes in the network. Each node, /, is 
assumed to have a generating utility, i.e., a set of generators, associated with it and an external demand. Clearly, the 
upper capacity of the generating utility, , at node / can be set to zero to model a demand center. For a generating 
node, the external load, d* t , can be set to zero. For transmittal or transshipment nodes, we set both generation and 
demand to zero. 

The transmission cost and capacity in the United States are posted on the World Wide Web as required by the 
Federal Energy Regulatory Commission. We assume that our system can access the transmission costs and capacities 
for the different regions of interest. 



Solution Approach 



The computer implemented tool described herein uses the probabilistic distributions of the electric load and of the 
cost function of each utility to build probabilistic distributions for the spot-market prices and for the electric trades in 
the system. Due to the use of discrete random variables, i.e., finite S t and K t one can consider all combinations of s 
and k at each hour. The problem of cost minimization while satisfying demand and profit margins is then solved for 
each one of these combinations to produce a spot-market price and to compute the trades. The results can be bucketed 
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example, R W. Cottle, J. S. Pang, and R. e. Stone, The Linear Complementary Problem, Academic Press, 1992, and 
K. G. Murty, Linear Complementarity, Linear and Nonlinear Programming, Helderman Verlag, 1988. 

Figure 2 shows a flow chart of the previous process applied to multiple utilities, 1 multiple scenarios, 1 S, 

and multiple periods, 1 T The process starts in function block 201 in which the input data is passed to this procedure. 
The input data are the historical prices, Wf t , Jf t , and the corresponding amounts of power generated, . 

These data are usually available from the power pool and assumed to be stored on the hard disk of our computer' In 
many cases, such data are available through the World Wide Web (WWW). At this stage, the counter / is set to one. 
The process then moves to function block 202 in which the counter s is set to one. Function block 203 sets the value 
of / to one. The process then proceeds to function block 204 in which the quadratic program of equation (1 ) is constructed 
and solved. We do not discuss the solution approach here. Most optimization packages provide specialized routines 
for handling such problems. We suggest using the Optimization Subroutine Library of I BM which provides a specialized 
function, EKKQSLV, for this operation. As a result, we obtain the values of af t , bf f and cf v which are used to approx- 
imate the price function, i.e., production plus profit margin, for any generation requirement. In function block 205, the 
value of t is incremented by one. Decision block 206 compares / with T. If t<T, the execution is directed to function 
block 204; otherwise, the execution proceeds to function block 207. In function block 207 the value of the scenario 
index s is increased by one. Decision block 208 checks if s is greater than S. If the answer is no, then the execution 
is directed back to decision block 203. Otherwise, the execution continues to function block 209. In function block 209, 
the value of / is incremented by one. Decision block 210 directs the flow to function block 202 if / is less than or equal 
to n\ otherwise, the process terminates. As a result, we obtain a set of quadratic functions that can be used to approx- 
imate the price functions under different scenarios and during different time periods. 

The Explicit Case 

The second case we discuss is that in whichlhe price function is known explicitly That is, for each generating unit 
owned by utility /; we know the operating range, the cost function within that range, the start-up cost, and any other 
characteristics relevant to a generator. We may also know the status of each generator. That is, we know if a generator 
is scheduled for maintenance and if it is unavailable due to mechanical problems. Furthermore, one may have historical 
information regarding the probability of failure of each unit. All of the previous information implicitly represents the cost 
structure of a utility That is, we do not have an explicit representation of If but we can compute the value of /? given 
the load and conditions of the generators. Note that in the current electric -power industry in the USA, the previous data 
are known to the general public. It is available on the World Wide Web and can be obtained in the form of a database 
stored on a floppy disk. However, due to competition, this information will not be available in the near future in which 
case we need to rely on the implicit representation mentioned above. 

To build a set of functions, f } that approximate the cost structure of a utility under different circumstances, we sort 
the generators in the increasing order of their average generating cost. This is done by adding the start-up cost of a 
unit to its operating cost at the maximum operating level and dividing by its maximum operating level. To clarify the 
previous approximation, the cost function of Figure 3A is used. To approximate this function, the start-up cost is added 
to the value of the function at the maximum operating level, G. Then, the average cost of this unit is computed by 
dividing the new cost at G by G. As a result, a generator has a linear cost function that passes through the origin and 
its operating range is anywhere between 0 and G (see Figure 3B). The generating units are sorted in the increasing 
order of their average costs which results in a convex piecewise linear function for operating the system of a utility. 
Figure 4 shows an example of such a function. The first, i.e., cheapest, generator is used to cover the load between 
0 and its 1 maximum capacity, G-,. The second most economic generator is used to meet the electric load between the 
additional load up to the capacity of the second generator. That is, it covers the load between G^ and G^+G 2 . The 
process is repeated until all the generating units of a utility are incorporated into the system. Note that we can build 
different functions to reflect the different operating conditions of a utility. For example, if some generators are not 
available, then we have a different function. Another example is the case in which a generator uses gas as its primary 
fuel. In this case, the cost function can be affected by gas prices in the market. Hence, we can create different functions 
(and different approximations) to model different market conditions. 

Note that the resulting function is an approximation of the cost function of a utility and not its price function. To 
obtain the price function from the cost function we need to estimate the profit margin of a utility The goal of any 
deregulated market is to force the competing entities to reduce their cost. The minimum cost at which a utility can sell 
its production, without losing money, is the production cost of its electricity However, this pricing scheme does not 
generate any profit to the utility. To avoid the problem of zero profit, each competitor usually sets a profit margin that 
is acceptable for him or her. This profit margin may be different for the different periods, t, of the planning horizon and 
under the various scenarios. For example, a profit margin of 1% may be acceptable in periods in which the market 
demand is low while a higher margin might be desirable if the demand is high. The margin can also change with the 
scenarios For example, if fuel prices are high, one may not be able to enforce high profit margins. The price function 
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the input parameters, data available about utilities and their generators, are not necessarily known in advance. To 
avoid this difficulty, we describe the input parameters in general terms so that our tool can handle any electric power 
market. The input to our system consists of (1 ) the generating capacity and price functions, ff t , of utilities competing 
in the market, (2) the load forecast in the regions of interest, and (3) the transmission costs, f , charged by the ISO 
or other equivalent body to transmit electricity from node / to node / at time period t and the transmission capacity, 
u^. t . Note that all of our notations contain the superscript s. This is the scenario index, a concept that is described 
in ifurther detail in the above referenced US patent application and which shall be clarified through the following dis- 
cussion. Here is a detailed description of the input parameters. 

1. The maximum generating capacity, G? r of each utility, /, competing within the market. Note that this capacity 
is not necessarily the total generating capacity of all generators owned by a utility. For example, some units may be 
unavailable due to maintenance during certain days. Also, depending on the season and availability of fuel, the capacity 
may change. The scenario index, s, allows a certain utility to have different possible capacities. The capacity could 
also change depending on rover the planning horizon, usually one week or 168 periods. For example, the capacity 
over the weekends is lower than regular days. 

2. The price functions, 1f t at which utility / is willing to sell its power at time t under sbenario s The value 
returned by f f * is assumed to include the operating cost plus the profit margin of a utility. Note that /? f may be given 
(a) implicitly or (b) explicitly. 

The Implicit Case 

We may not know the exact cost structure of a utility but may be able to approximate its operating cost by a 
quadratic function. For instance, we may observe that a utility /, priced its power at time t and under a certain load at 
Xf $/MWH. By observing the price offered by a utility for its power over time, one may construct an approximation of 
the cost structure of this utility. The simplest way to do so is to collect the different data points of price versus load and 
fit a quadratic curve using least squares. Our experience indicates that the quadratic approximation is a good one for 
most utilities. Of course, if one observes a certain profile for the cost structure of a utility, one may consider a different 
function to be fit. 

The problem of "guessing" the price function of a utility is not a simple one in a competitive environment. For 
example, one may observe a cost of $15 per MWH for a certain utility under a 12,000 MWH load on one day while 
observing an operating cost of $18 per MWH for the same utility under a similar load in another day The reason for 
this discrepancy might be a generating unit that is going through maintenance, a failure of a generator in that utility, a 
large error in the load forecast which resulted in shortage of generation, or a shortage of certain fuel which may have 
forced the utility to use a more expensive fuel. Regardless of the reason, building the price function ff of a utility / may 
be impossible if one has to use a single value for the cost under a given load. To handle this difficulty, our model permits 
the use of more than one price function, f* , at any time period t For a given time period f, we allow the price function 
to be any of the functions f] { . ff t where S, is the number of different functions available at time t. That is, we model the 
cost function of a utility / at each time period r by a set of cost functions. The index s is read as scenario. One hopes 
that using a set of functions reflects accurately the uncertainty an observer may have regarding the operating cost of 
a utility. For example, one scenario may represent a normal operating environment, another might be the cost observed 
when natural gas prices are high, a third may be the price function when there is a heat wave in the area, and so on. 
One can use as many scenarios as needed so that the uncertainty in / is approximated reasonably 

Given a set of observations, X J ; s r /=1 . ■■■Jf r . of the prices of utility /at time t under certain conditions s, we find the 
pricing profile that best-fits this utility using a least square approach. Here is a pseudocode of this process. 

Initialization. Get a set of prices, )Jf t , offered by utility /at time r under scenarios. Let ^ be the amount of electricity 
generated by the electric utility for these price offers. 

Least Squares. The approximating function must be quadratic, convex, positive, and increasing. To find the function 
that best-fits our data while satisfying the previous constraints, we solve the following quadratic program: 




subject to 




Kc ±0, 



The previous optimization model can be solved efficiently using any quadratic programming approach. See, for 
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Table 1 (continued) 
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is The first column in the above table represents the spot market price. The second column contains the amount of power 
that may be traded. Note that for the same spot price, one may buy or sell different amounts of electric power. The 
likelihood of a transaction is given in the last column. 

It will be appreciated that it is desirable to use all available historical, competitive, and weather information to 
produce the most accurate output, although the data that is actually available for different utilities may differ from one 

20 region to another. For example, some utilities are required, by the power pool, to post the status of their generating 
units by the end of the week. These requirements differ significantly between power pools. Also, data that are currently 
public, such as hourly loads, may not be available in the near future. Accordingly, the tool can be operated in a flexible 
manner in accordance with the amount of information that is in fact available. 

A preferred embodiment of the invention will now be described in detail by way of example only with reference to 

2S the following drawings: 



Figures 1 A and 1 B are block diagrams illustrating the relation of the power trading and forecasting tool to its inputs 
and outputs; 

Figure 2 is a flow diagram illustrating the process of approximating the price function of a utility by a quadratic 
30 function; 

Figures 3A and 3B are graphs demonstrating how the cost function of a generator can be approximated using a 
linear function; 

Figure 4 is a graph showing an approximate cost function of a utility constructed by using the linear approximations 
of Figure 3; 

35 Figure 5 is a flow diagram of the simulation used in the process; 

Figure 6 is a flow diagram of a process that generates a random sample; and 

Figure 7 is a flow diagram showing the bucketing process used to create a probabilistic distribution using the output 
of simulation runs. 

40 This application is related in subject matter to U.S. Patent Application Serial No. 08/869,561 (IBM Docket 

YO997-032) filed June 5, 1997, by Samer Takriti for "RISK MANAGEMENT SYSTEM FOR ELECTRIC UTILITIES". 
This earlier application is incorporated herein by reference, and a copy of this earlier application is available on the file 
of the present application. 

Referring now to the drawings, and more particularly to Figures 1 A and 1 B : there are shown in block diagram form 
45 the relationship of the power trading and forecasting tool and its inputs and outputs. The tool can have different sets 
of input parameter and two are described here by way of example. The first example is shown in Figure 1 A. In this 
case, we are given the historical data of selling prices and amounts of power sold by each utility. The tool 11 receives 
inputs from a plurality of utilities ^2 A to 12„. These inputs include for each utility prices at which it sold power in the 
past and the corresponding amounts of power sold. In addition, the tool 11 receives as an input the transmission costs 
50 and the maximum transmission capacity for different lines in the transmission network 13. The output of tool 11 is the 
probabilistic distribution for trades and spot prices 14. 

The second set of input parameters is illustrated in Figure 1 B. They include the operating cost curves of the different 
generators 12j to 12„, estimates of the marginal profits, and electric load forecasts. The output is similar to that of 
Figure 1 A. The inputs and outputs of the tool 11 are described in more detail below. 

55 

System's Inputs 



Since the details of the deregulated market may vary from country to country, and so cannot be completely defined, 
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probability for each range; 

and generating a histogram of the probabilities of each range as an approximation of the probability distribution 
of the marginal cost at each delivery point at a given time. 

5 The invention further provides a system for use in an electricity network to manage power supply by forecasting 

spot-market prices of electric power and trading transactions at different delivery points, said system comprising: 

means for inputting maximum generating capacities of each utility competing within the market, price functions at 
which a utility is willing to sell its power at a given time period, a weather forecast, and a transmission cost and 
10 transmission capacity in different power lines of the electrical system; 

means for using probabilistic distributions of electric load and the price function of each utility and transmission 
cost and transmission capacity between different utilities to build probabilistic distributions for the spot-market 
prices and for electric trades in the system; and 

means for outputting probabilistic distributions for spot prices and trades at a plurality of time periods of a planning 
is horizon and for a plurality of delivery points to assist in said supply management. 

Viewed from another aspect the invention provides a computer power trading and forecasting process for fore- 
casting spot -market prices of electric power and trading transactions at different delivery points comprising the steps of: 

20 inputting maximum generating capacities of each utility competing within the market, price functions at which a 

utility is willing to sell its power at a given time period, a weather forecast, and a transmission cost transmission 
capacity in different power lines of the electrical system; 

using probabilistic distributions of electric load and the price function of each utility and transmission cost and 
transmission capacity between different utilities to build probabilistic distributions for the spot^market prices and 
25 for electric trades in the system; and 

outputting probabilistic distributions for spot prices and trades at a plurality of time periods of a planning horizon 
and for a plurality of delivery points to enable a power broker to effectively manage risk in buying and selling electric 
power with different suppliers and consumers. 

30 in a preferred embodiment, the probabilistic distributions of the transmission costs and capacities of power lines 

are created from historical data, and the probabilistic distributions for the spot-market prices and for electric trades in 
the system are created by sampling from electric load, price functions and cost distributions to create a sequence of 
deterministic models. Preferably, the process further comprises the steps of minimizing the cost of operating an electric 
utility system under each sampled set of price functions, transmission costs and electric loads; and storing a resulting 

55 output. The sampling from electric load, transmission costs and cost distributions is performed using a random number 
generator. 

Thus the computer implemented process described herein forecasts the spot-market prices and the trading trans- 
actions at different delivery points, providing the decision maker with probabilistic distributions for spot prices and 
trading to allow efficient risk management. The user of the process can be (1 ) a utility or an independent power supplier 

40 who is interested in using the trading and spot market information to decide which generating units to commit, how to 
price power, and the level of generation needed; or (2) a power broker who is interested in buying and selling power 
(but not producing) in order to maximize revenue. The computer implemented process uses the weather forecast in 
addition fo other market information and provides probabilistic distributions for the power traded and spot-market prices, 
rather than a single estimate for the spot price and the power traded. An example of the output of the computer imple- 

^5 mented process is shown in Table 1 below; the information therein can be used for hedging against uncertainty thereby 
allowing risk to be acceptably contained. 
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ZXi^J^S^^ "nes as well as other operational factors. To maintain a reliable 

tosubmithH.^ n ° UnCeS : , T ^ ■ ,OUrhours,na ^ ancetheloadfore ^t on the system and asks interested suppliers 
to submit b ds. ..e.. amount of power to be delivered at each point and pricing information. The load forecast is oiven 
for each deNve-y point within the physica. system. Bids are submitted to the Power Exchange orTx Shen K 
an aucon to de.erm.ne whfch suppliers to buy power from. For each delivery point, the PX cheeses.^ 

yZ^mTNo e e tharrr d ^ *" " " ^ ^ iS ' epea,ed ^ the 

system is met. Note that such a process ensures that the customers' demand for electricitv is satisfied at * minim*. 

St ^sr* °r supp,ier f ch ° sen * *• ^ *• s^^s^t^ ic 

point The PX performs this process for each delivery point in the system and at each hour of the day 

After agreeing on the production schedule, the ISO studies the schedule to make sure that it is feasible i e it 
does not exceed the capacrty of the transmission lines. Transmission charges are a function oHhe c^nqes'on of the 
power lines. These charges are calcu,ated by ISO and bii.ed to the different parties ^S^TSS^J^ 

™2T» * 3 W ' ,h 3 hi9h PMCe 6nd UP n0t S6llin 9 anv of their P"***»- On the other htnd seZ 
Tolt 2, PnCe T ^ 'T 18 en ° U9h reV6nUe f ° r 3 9en0ratOr Given ,hat no °™ knows in advancement 
£fZ , -"T T e,ectric -P° wer mar ket will become more uncertain and risky The hope is 

that compet,t,on w I force utilities to minimize their production costs whk:h will, in turn, benefit consumers and Ihe 

TZZZ 9 Th ' S m0de ' alS ° enCOUra9eS individualsl ° -est money in the electric-power marke Fo nstance 

oowei ^Sn7in S tn C r * "tT "* ***** ^ l ° C ° n9eSti0n ' an inves <° r mav ^ a " a ItematS 

power l,ne ,n the hope of benefitting from the high-traffic in that region. The result, of course, is a reduction in the 
transmission cost ,n that region and an increase in the reliability of the system reaucnon ,n the 

nuJirlr 6 T Pr , e ,^ ed ' hat USA iS 001 a ' 0ne in introducin 9 deregulation. For example, in the UK the formerly 
pub be electncty utility (a monopoly) has been split at privatisation into competing generators and multiple distrZos 
Again, wh st such deregulation is expected to lead to lower prices, it signLn.ly'adds to 

in the electricity network. In particular, it will be recognised that in the face of deregulation, the propeTZnagZToi 
TZ^Z P Tr ?l 3CCUrate f ° reCaS,ing " ° rdSr to aVOid ,6Chnical and/or commeLnai,rreTorme 

bv fa^Sf th t 6 inV ! n ! i0n Pr ° VideS 3 C ° mpU,er pr0C6SS for use in an elec,ricit y ne,wo * to manage power supply 
cZ^Z Zso! PTOeS °' 6leCtriC P ° Wer ,radin9 tranSaC,i ° nS * different deliver * Points 'said meS 

inputting maximum generating capacities of each utility competing within the market, price functions at which a 

caoar'^/'' h P^"^ a t a given time period, a weather forecast, and a transmiss^cfl^os^a^d transmission 

capacity in different power lines of the electrical system- transmission 

using probabilistic distributions of electric load and the price function of each utility and transmission cost and 

!=»cr2r uti,ities to build probabi,is,ic ~ s f - - ™° s s 

«SKSS^ of <~ - a *-» — 

Such a forecasting process provides the technical and commercial information necessary for producinq or distrib 

ri!«!t h' C »^ er 0 veasa '' s ' ac ' or y °P era '' orl of the electricity network, and the provi^on of^o^p^ut^o^a^bnlstk; 

distnbuttons in particular helps in effective risk managment. P proDabllls,lc 

It is preferred that the process further comprises the step of defining different scenarios wherein the maximum 
generating capacity o, each uti.ity changes under said different scenarios and also furthe cc^pTse h sip liuZ 
the input weather forecast to create a probabilistic distribution of load forecasts for each utility of interest 9 
QuadraticTnL P rwr h C,i0nS * f ° r 6Xample by ^PP^imating a price Junction of a utility as a 

ESS^'5S£^^ ^ inCreaSin9 39ainSt ^ ° r ~> be g.en explicitly as a 

sJL^lZTJT^T probabi,istic distribution for the spot-market prices and for e.ectrb trades in the 
terirm^ I S r P " n9 ,r ° m 6leCtriC '° ad ' priCS ,Unctions and cost distributions to create a sequence of de 

S Zi , ' 6 C ° St ° f ° Perating 3n 6,eC,riC Uti,ity svstem bein 9 noised under each samplil se of price 
functions, transmission costs and electric loads. P 

buJ«Z er ^ Pref6rred embodiment ' the Probabilistic distributions for spot prices and trades are determined by a 
bucketing process comprising the steps of: oioiminaa oy a 

creating ranges of spot prices; 

counting the number of spot prices that fall within each range and normalizing the counted numbers to obtain a 
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Description 

The present invention generally relates to a computer process and system for use in an electricity network to 
manage power supply by forecasting spot -market prices of electric power and trading transactions at different delivery 

s points based on various technical and commercial criteria. 

Electricity is an essential part of our lives. If we took a moment to think of the services that would not exist if 
electricity were not available, we would be amazed. Almost every aspect of our modem lives involves electric power, 
from light bulbs and television sets to hospitals and automobile industries. Although we are used to having power 
whenever we need it, the processes and systems involved in delivering electricity require careful planning and sophis- 

10 ticated mathematical and forecasting models, in order to allow for various technical, environmental, and commercial 
considerations. Examples of technical considerations include the operating characteristic of power stations (which may 
not be able to respond instantly to demand changes), and capacity limitations of the transmission network; an example 
of an environmental consideration is the weather, which has a strong impact on consumer demand for electric power; 
and an example of a commercial consideration is the action of a competitor (eg the amount and price of electricity 

is which it is prepared to supply). 

Lately, due to rising costs of energy and discrepancies in its price among different regions of the country, the legal 
framework surrounding the electric-power industry in the USA has changed (a process known as deregulation). This 
change has opened the door for utilities to compete with each other and against independent suppliers regardless of 
their geographic location. Although this change will benefit the consumer, utilities are going to face a highly unpredict- 

20 able market and wilt need to make tough decisions regarding power generation and delivery. The current picture of a 
single utility controlling the market in a specific region will soon disappear Instead, there will be power producers who 
sell their product to a power pool; and power suppliers who will buy power from this pool and in turn sell it to their 
customers. Although the full picture of the power industry after deregulation is not yet known, it is clear that utilities 
need to prepare themselves for an open market in Which buying and selling power are to be considered when scheduling 

25 their generating units. The main reason behind such deregulation is to reduce the high price of electric energy. Initial 
steps towards deregulation were taken in 1978 with the passage in the USA of the Public Utilities Regulatory Policy 
Act. This act encouraged non-utility generation and required utilities to buy power from independent generators. The 
Energy Policy Act of 1992 took deregulation a step further by mandating open access to the transmission system for 
wholesalers. 

30 Currently, electricity is sold as a service that is delivered to specified points. For example, one typically expects to 

receive electric power via a meter outside the house. We pay for this service regardless of its producer or which power 
lines it followed. That is, an electricity bill indicates the total usage of electricity in kilowatt hours (KWH) and the service 
price per KWH without incorporating any other details into the pricing scheme. Deregulation in the USA is changing 
this picture by unbundling the electric power into generation and transmission. In the future, one will pay a production 

35 cost and a transmission fee. There will be several power suppliers from whom electric power may be purchased. 
Suppliers may have different pricing mechanisms. For instance, there might be a discount for using power off-peak 
periods or for signing a long-term contract with the supplier. Power producers will compete with each other to minimize 
their costs so that they can sell their product to more customers and maximize their profit. 

The hope is that deregulation will result in cheaper prices and play a part in improving the economy by encouraging 

40 investment in electric utilities. The size of the electric industry is expected to grow after deregulation as was the case 
with the telecommunications industry. The telecommunications industry's revenue shot up from $81 billion to $170 
billion within ten years of deregulation. ; 

A pitfall of deregulation is that the load on a utility system is becoming increasingly unpredictable. The reason is 
that trading transactions can change the load pattern significantly For example, some utilities may sell more than 30% 

45 of their power generation to other utilities on certain days. Demand and supply in the market are functions of volatile 
electricity prices which in turn depend on highly unpredictable elements such as regional weather conditions and fuel 
prices. 

There are two components for any competitive environment: physical and commercial. In the case of electricity, 
the physical component already exists in the form of transmission lines and all consumers have physical access to the 
50 market. The commercial component consists of the ability of buyers and suppliers to execute trades over the physical 
network. One should also ensure that anyone who is interested in joining this market has nondiscriminatory access to 
the physical network. In this environment, suppliers and consumers have the freedom to contract with each other under 
negotiated terms and conditions. To ensure the delivery of the power traded, they need a mechanism to ensure the 
physical delivery of the electricity Also, they may need to rely on a spot market to buy or sell additional electricity. 




Although the new environment in the USA encourages competition, power transmission will remain regulated there 
for the time being in order to maintain a reliable system. The transmission lines in each state or region will be controlled 
by an independent entity called an Independent System Operator or ISO. One of the ISO responsibilities is to settle 
financially with the parties involved in transmitting electric power. The transmission cost depends on the proximity of 
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